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Abstract 

We present an analytical calculation of the next-to-next-to- leading order correc- 
tions to the partial decay width H ^ gg for a Higgs boson in the intermediate mass 
range. We apply an asymptotic expansion for Mjj <C 2Mt and compute three terms 
in the expansion. The leading term confirms the results present in the literature. It 
is argued that our result is equivalent to an exact calculation up to Mh ~ Mt. For 
a Higgs boson mass of 120 GeV the power-suppressed terms lead to corrections of 
about 9% in the next-to-next-to-leading order coefficient. 

PACS numbers: 14.80.Bn 12.38.-t 



1 Introduction 



Up to date the Standard Model (SM) Higgs boson has evaded its direct detection. Elec- 
troweak precision data collected at the CERN Large Electron- Positron Collider (LEP), 
at the Stanford Linear Collider (SLC) and at the Tevatron at Fermilab predict a light 
SM Higgs boson below approximately 200 GeV. On the other hand, the direct search at 
LEP2 has excluded Higgs bosons below 114 GeV which leaves a relatively narrow window 
for the mass. A Higgs boson in this mass range, often also referred to as intermediate- 
mass Higgs boson, predominantly decays into bottom quarks and W bosons depending 
whether the latter decay is kinematically allowed. A further decay channel which is of 
phenomenological interest is the one into gluons. At lowest order this process is mediated 
by a top quark loop [1,2]. The next-to-leading order (NLO) corrections are quite large 
and amount to roughly 70% [3-8]. About ten years ago the NNLO corrections have been 
computed in the limit of an infinitely heavy top quark mass [9,10]. Since the corrections 
amount to approximately 20% they increased the reliability of the perturbative expansion 
of the decay rate. This has been further substantiated by the recent evaluation of the 
NNNLO corrections [11], again in the infinite top quark mass limit, which provides a 
contribution of about 2%. 

The calculations to NNLO and NNNLO are both based on an effective theory which results 
from the SM after integrating out the top quark. The effective Lagrangian is constructed 
from operators with dimension five which is sufficient to obtain the leading term in the 
heavy-top expansion. The inclusion of higher- dimensional operators would lead to power- 
suppressed terms. However, the renormalization of the corresponding effective Lagrangian 
is quite tedious. In this paper we follow a different approach and consider the Higgs-boson 
propagator in the full SM. Its imaginary part leads — with the help of the optical theorem 
— directly to the decay rate. The disadvantage of this method is that at NNLO five- 
loop corrections to the Higgs two-point function have to be considered since the LO 
result already requires a three-loop diagram where the coupling of the Higgs boson to 
the gluons is mediated by two top quark triangles. They are connected by two gluon 
lines. The advantage of our method relies on the straightforward automatization using 
state-of-the-art techniques and program packages as we will show in Section [2l 

An important feature of the Higgs boson decay into gluons is its affinity to the production 
mechanism in the gluon fusion channel. For the latter the NNLO corrections are again only 
known within the framework of the effective theory [12-14]. The calculation performed 
in the present paper gives a first hint about the potential size of the terms which are 
suppressed by the top quark mass. 

Let us finally mention that the corresponding NNLO corrections to the decay of the 
Higgs boson into two photons have been obtained in Ref. [15]. Note that in this case a 
simple expansion of the Higgs-photon-photon vertex diagrams is sufficient to obtain power- 
suppressed terms the heavy Mf-limit. Thus only three- instead of five- loop diagrams have 
to be considered. 



2 




Figure 1: Sample diagrams contributing to Hh- The curly lines represent gluons and the 
solid lines stand for light quarks or top quarks. The Higgs boson (external dashed lines) 
only couples to top quarks. 



The paper is organized as follows: In Section [2] we present some details of our calculation 
and Section [3] contains our analytical results and the discussion about the numerical 
implications. We conclude the paper with Section |H 

2 The calculation 

As already mentioned in the Introduction, for the computation of the decay rate r{H 
gg) we consider the Higgs boson self energy, 11// (g^), and apply the optical theorem which 
in our case reads 

TiH^gg) = ^Im [n//(Mi + zO)] . (1) 

Since we do not consider the (exclusive) decay into light quarks only those diagrams 
contribute to Uh where the Higgs boson couples to the top quark. Thus the LO result is 
obtained from the three-loop diagram shown in Fig.[T]^a) (and the one with crossed gluons). 
The NLO and NNLO results are obtained by dressing this diagram with additional gluons 
and (light and heavy) quark loops. Sample diagrams can be found in Fig. [T](b)-(e). It is 
clear that beyond LO next to gluons also light quarks can be produced. In particular, at 
NNLO also a final state qqq'q' with light quark flavours q and q' is possible. 

Our approach for the computation of r{H gg) requires the evaluation of five- loop dia- 
grams in order to obtain the NNLO corrections. With the currently available techniques 
their exact evaluation is out of reach. However, for a Higgs boson in the intermediate 
mass range it is promising to consider an expansion for Mfj <C which shows a rapid 
convergence at LO and NLO as we demonstrate in Section [31 In this limit we can apply 
the so-called hard- mass procedure (see, e.g., Ref. [16]) to the propagator diagrams which 
then factorize into products of one-, two- and three-loop integrals. A graphical example 
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Figure 2: Sample five-loop diagram contribution to Hh where the hard-mass procedure 
is applied in graphical form. Only the subgraphs contributing to the imaginary part are 
shown. 



is shown in Fig. [2] where — according to the rules of the hard-mass procedure — five 
subdiagrams are identified. They have to be expanded in the small quantities and then 
to be inserted in the co-subgraph which (in this example) only consists of gluons. In 
Fig. [2] they are sandwiched between the symbols 'V . The hard-mass procedure has been 
automated in the program exp [17, 18] which we use for our calculation. 

At LO we have two contributing Feynman diagrams, at NLO there are already 71 and 
at NNLO 2649. All Feynman diagrams are generated with QGRAF [19] and the various 
topologies are identified with the help of q2e [17,18]. Afterwards exp [17,18] applies the 
hard-mass expansion and produces output which can be read by the FORM [20] packages 
MATAD [21] and MINCER [22]. MATAD can compute vacuum integrals up to three-loop order 
where in our case the mass scale is given by the top quark. On the other hand, MINCER 
can handle massless propagator-type diagrams where for our application the external 
momentum is given by the Higgs boson mass. 
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According to the optical theorem we have to compute the imaginary part of Uh which 
requires the evaluation of the pole part of the three-, four- and five- loop diagrams. Since 
MATAD and MINCER are constructed such that at most the finite parts of three-loop diagrams 
can be computed slight modifications are necessary which allow the correct evaluation of 
the 1/e pole part of the five-loop diagrams. 

For some of the diagrams it turned out to be crucial to use the parallel versions of FORM 
— ParFORM [23] and TFORM [24] — in order to obtain results in a moderate amount of 
time. With our setup, which is not tuned for a parallel computation, we obtain speedup 
factors between four and five using eight processors. 



3 Results and discussion 



3.1 LO and NLO result 



The method described in the previous Section can be tested at LO and NLO where the 
exact results are known. Let us for this purpose introduce the notation 



TiH-.gg) = Lq 



TT 



pLO _|_ Q^sIA^) A^-pNLO 



+ 



7r 



Ar 



NNLO 
99 



+ 



,(2) 



with To = GFM^/(367rv^). F^o is given by [1,2] 



gg 



2r 



1 ^ arcsin^ a/t 

226 . 
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(3) 



with r = Mfj/{4M^) where Mf is the on-shell quark mass. In the second line of Eq. 
we have performed an expansion up to order which we could confirm by applying the 
hard- mass procedure to the three- loop diagrams (cf. Fig [I](a)). The ellipses indicate 
higher order terms in r. 

The exact NLO result is only known in numerical form [4-8]. Thus we concentrate on the 
expansion for Mjj -C 2Mf. The first three terms have been computed in Ref. [25] which 
we could confirm. Furthermore, we have added two more terms and arrive at 
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where 
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with Lh = \n{n'^/Mjj) and Lt = \n{ij^/m'f). In Eq. (jl]) we have f = M|^/(4m^) where 
nit = f^tif^) is the MS top quark mass. For convenience we have kept the label referring 
to closed light-quark loops which takes the numerical value ni = 5. Furthermore, we have 
asifi) = a'f\n), i.e., the top quark has been decoupled from the running of the strong 
coupling (see, e.g., Ref. [10]). 

Transforming the top quark mass from the MS to the on-shell scheme [26] leads to the 
following result which is evaluated in numerical form 

^pNLO ^ 17.9167 + 3.83333Lj^ + (9.574+ 1.789LH)r 
+ (5.571 + 0.939Lh) + (3.533 + 0.550Li^) 
+ (2.395 + 0.350Lh) + (1.708 + 0.237^^^) + . . . . (6) 

In Fig. [3] is shown as a function of r (solid line) and compared to the approximations 
where successively higher orders are included (dashed lines). The analog plots for the 
NLO results Ar^j"° and Ar^^^° are shown in Fig. H] where /i^ = has been chosen. 
From the behaviour of the approximations we can deduce that in the region, where the 
curves including different orders in the M^-expansion are on top of each other, they also 
coincide with the exact result. Having this in mind we conclude for the NLO corrections 
that in case of the pole mass the approximations including the (r^) terms coincide with 
the exact result up to r ^ 0.60 (0.35) which corresponds to Mh ~ 270 (200) GeV. It is 
worth to mention that the top quark mass-suppressed corrections become smaller if the 
MS mass is used for the parameterization. In this case the curves including the f ^ and 
terms, respectively, are basically on top of each other — almost up to f = 1. Even 
the curve including corrections up to order shows a deviation of less then 6% for f = 1 
and provides a perfect approximation up to f = 0.5 (i.e. Mh ~ 250 GeV). Thus, it can 
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Figure 3: as a function of r where successively higher orders are included (dashed 
lines). The exact result is shown as a solid line. 
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Figure 4: AT^^^ (left) and AF^J"*-* (right) as a function of f and r, respectively, where 
successively higher orders are included. For the renormalization scale /s^ — Mjj has been 
chosen. 



be assumed that in the phenomenologically relevant region for Mu the results including 
the first three expansion terms in r represent an excellent approximation for all practical 
purposes. Let us also mention that for Mh = 120 GcV (corresponding to r ^ 0.1) the 
leading term in the 1/Mt expansion approximates the exact result with an accuracy of 
about 5% — both for the on-shell and MS top quark mass. These considerations provide 
a strong motivation to compute the first three expansion terms at NNLO. 

Very often, in particular in the context of the Higgs boson production in the gluon fu- 
sion process, the complete LO result is factored out when considering the perturbative 
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Figure 5: ^^^^^ (left) and ^^^^^ (right) as a function of f and r, respectively, where 
successively higher orders are included. For the renormalization scale fi"^ = M'jj has been 
chosen. 



expansion. In this case we define 



T{H^gg) = T,[^^] T 



TX 



2 



99 



a.(/i) . NLO , ( ^^^(aO \ ^ epNNLO , 

^ 99 I ^ J " 99 "T ■ • • 
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The coefficients 5r?]"° — both in the MS and on-shell scheme — can easily be obtained 



99 



from the results for AT^^^. We refrain from listing them explicitly. 

In Fig. [5] 5r^j"° and ^r^^^*-* are shown as a function of f and r, respectively, where 
successively higher orders are included. One observes a dramatic improvement in the 
convergence of the l/Mf expansion for r <^ 0.3 where already the leading Mf-term is 
practically indistinguishable from the exact result. 



3.2 H ^ gg to NNLO 

Applying the methods described in Section [2] to the five-loop diagrams shown in Fig. [1] 
leads to the NNLO corrections to the decay rate T{H gg). We were able to compute 
three expansion terms in r which can be cast in the form 
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where the MS top quark mass has been used. Transforming the result to the on-shell 
scheme and evaluating it in numerical form leads to 

^pNNLO ^ 156.808 + 102.146Lh + 11.021L^ + 5.708Lt + 
(109.365 + 52.662Lh + 5.143L^ + 4.645Lt) t + 
(74.434 + 29.920Ljf + 2.699L| + 3.297Lt) + . . . , (10) 

with Lt = ln{fi'^/M^). We want to mention that the result for h^'^^" has already been 
obtained in Ref. [9, 10], the remaining terms are new. 

Let us for completeness mention that the leading term of the large top quark mass ex- 
pansion is also known at NNNLO [11]. It is given by 

^pNNNLO ^ 467.684 + 1215.302 Lh + 394.514 + 28. 164^3^ 

+ 21.882 LhLt + 122.441 Lt + 10.941 L| + . . . . (11) 

In Fig. [6] we show the results of Eqs. (IE]) and ( fTOl) . For small values of r one observes a 
logarithmic divergence which even occurs for fi = Mh- In the phenomenologically relevant 



9 





Figure 6: AF^j^^^ (left) and AF^j^^*-* (right) as a function of f and r, respectively, where 
successively higher orders are included. For the renormalization scale /z^ = Mfj has been 
chosen. 



range (r ^ 0.1) the logarithm In (M^/M^^) is numerically small and we observe the same 
pattern as at LO and NLO which is a strong indication that our result including the 
terms provides an excellent approximation to the exact result up to Mh ~ Mt {r ^ 0.25). 
For Mh = 120 GeV (r ^ 0.1) we observe a correction of about 9% originating from the 
power-suppressed terms. This induces an effect of approximately 1% on the decay rate 

Fig. [7] shows the result for (JF^^'"^ and (JF^'"^'^'-'. Similarly to the NLO result one observes 
a reduction of the size of the mass corrections. This is particularly pronounced in the 
case of the MS scheme where the curves involving terms up to order and lie on 
top of each other and coincide with the leading term up to f ^ 0.4 {Mh ~ 220 GeV). 
These considerations provide a strong motivation to consider the total decay rate where 
the complete LO result is extracted. Furthermore, due to the similarity to the production 
of Higgs bosons in the gluon fusion channel it can be expected that there the mass 
corrections are also at the few percent level — below the current uncertainties from the 
parton distribution functions [12-14,27]. 

Let us finally mention that the numerical effect of the power-suppressed terms computed 
in this paper are comparable to the leading Mj-term at NNNLO (cf. Eq. ( fTTi) ). 

At this point it is interesting to consider the dependence of r{H gg) on the renor- 
malization scale /i. In Fig. [S] the LO (dashed), NLO (dash-dotted), NNLO (solid) and 
NNNLO (dotted) result is shown for fi between 10 GeV and 1 TeV where for the LO 
curve the exact result has been used. The NLO and NNLO curve include terms of or- 
der and r^, respectively. As input for our numerical studies we use as{Mz) = 0.118, 
Mz = 91.1876 GeV [28], Mt = 170.9 GeV [29] and Mh = 120 GeV. The renormalization 
group equations are solved with the help of the program RunDec [30]. 

It is interesting to note that around /i = 40 GeV the NNLO curve shows a local maximum 



10 



220 

200 

2 180 - 

2 OX) 

'fe 160 - 
140 
120 

0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

T T 

Figure 7: (5r^j^'"° (left) and 5r^j^'"° (right) as a function of f and r, respectively, where 
successively higher orders are included. For the renormalization scale fi^ = M'jj has been 
chosen. 



where the decay rate is independent of /i. For /x = 50 GeV the NNLO corrections vanish 
(crossing between dash-dotted and solid line) which hints for a good convergence of per- 
turbation theory. This is further supported by the smallness of the NNNLO corrections 
(dotted) which show a very flat /i-dependence over the whole range considered in Fig. [81 

4 Conclusions 

In this paper the NNLO corrections to the gluonic decay width of the Higgs boson in 
the intermediate mass range is considered. With the help of an automated asymptotic 
expansion three terms in the large- M^-limit have been obtained. It is demonstrated that 
our result is equivalent to an exact calculation at least up to Mh = Mt. For Mh = 
120 GeV the NNLO term which changes by about 9% due to the new corrections induces 
a 1% correction to the total gluonic decay rate. The situation is different if the complete 
top quark mass dependence of the leading order result is factored out. In this case the 
power-suppressed corrections become smaller. This is in particular true for top quark 
masses in the MS scheme where for Higgs boson masses up to 220 GeV the full result 
is given by the leading term in the large Mj-expansion. This observation has possible 
implications for the Higgs boson production via gluon fusion where only the leading order 
term in Mn/Mt is available. In case the expansion for Mh -C Mt shows a similar behaviour 
the uncertainty induced by the unknown power-suppressed terms is negligible. 

Acknowledgements 

We would like the thank Robert Harlander for useful comments and carefully reading the 
manuscript. This work was supported by the BMBF through 05HT6VKA and by the 
DFG through SFB/TR 9. 




11 




10 

II (GeV) 



10 



Figure 8: Tgg/To = T{H gg)/To as a function of n with Fq = GFMfj/(367r^/2). The 
LO, NLO and NNLO results are shown as dashed, dash-dotted and sohd hues. The dotted 
curve also includes the leading top quark mass dependent term at NNNLO [11]. 
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